wmo®  BV  tell/.  4  04806 


THERMOELASTIC  DAMPING 


t  TECHNICAL  DOCUMENTARY  REPORT  NO.  ASD-TDR-62-1031 


DECEMBER  1962 


DIRECTORATE  OF  MATERIALS  AND  PROCESSES 
AERONAUnCAL  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


O  DC 

Pniect  No.  7351,  Ta*  No.  735106 


TISIA  D 

(Picpand  tmder  Contnct  Na  AF  33(657)-7453 
by  die  Unhcnity  of  MiniieMila,  Minnrindii,  Mboii^ 

L.  E.  Goodman,  C.  S.  Chang,  A.  R.  Robmsoo,  andion) 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than 
in  connection  with  a  definitely  related  Government  procurement  operation,  the  United  States  Gov¬ 
ernment  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever;  and  the  hurt  that  the  Gov¬ 
ernment  may  have  formulated,  furnished,  or  in  any  way  supplied  Uie  said  drawings,  specifications, 
or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufiicture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  requesters  may  obtain  copies  of  this  report  from  the  Armed  Services  Technical  Infor¬ 
mation  Agency,  (ASTIA),  Arlington  Hall  Station,  Arlington  12,  Virginia. 


This  report  has  been  released  to  the  Office  of  Technical  Services,  U.  S.  Department  of  Com¬ 
merce,  Washington  25,  D.C.,  for  sale  to  the  general  public. 


Copies  of  this  report  should  not  be  returned  to  the  Aeronautical  Systems  Division  unless  re¬ 
turn  is  required  by  security  considerations,  contractual  obligations,  or  notice  on  a  specific  document 


FOREWORD 


This  report  was  prepared  by  the  Department  of  Aeronautics 
and  Engineering  Mechanics  of  the  University  of  Minnesota  under 
USAP  Contract  No.  AP  33(657)“7^53»  This  contract  was  initiated 
under  Project  No.  7351*  "Metallic  Materials,"  Task  No.  735106, 
"Behavior  of  Metals."  The  work  was  administered  under  the 
direction  of  the  Metals  and  Ceramics  Laboratory,  Directorate 
of  Materials  and  Processes,  Deputy  for  Technology,  Aeronautical 
Systems  Division,  with  Mr.  D.  M.  Pomey,  Jr.  acting  as  project 
engineer. 

This  report  covers  work  performed  during  the  period 
June  15,  1959  to  September  1961. 

The  authors  are  Indebted  to  Mrs.  Gerald  Webers  who 
typed  the  manuscript  and  to  Mr.  George  Peng  who  prepared 
the  figure. 


ABSTRACT 


The  non-linear  coupled  field  eqtiatlons  of  thermoelastlclty  are 
herein  solved  for  the  quasl-statlc  behavior  of  a  solid  bounded  by 
two  parallel  planes.  The  mechanical  energy  converted  to  heat 
dxirlng  each  cycle  of  the  loading  process  is  explicitly  evaluated 
and  Its  spectral  variation  computed.  From  this  exact  solution  an 
approximate,  but  quite  general,  expression  for  the  thermoelastic 
energy  dissipated  In  solids  as  a  result  of  elastic  deformation  Is 
developed . 
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I.  INTRODUCTION 


Therooelastlc  damping  as  one  of  the  cauaes  of  energy  dissi¬ 
pation  In  engineering  materials  has  been  studied  by  a  nuwer  of 
authors.  Anong  them  Zener  (1,2),*  Lessen  (3,4),  and  Blot  (5) 
may  be  mentioned.  In  his  treatment  of  thermoelastlclty  as  a  part 
of  general  theory  of  Irreversible  thermodynamics.  Blot  works  out 
an  example  of  thermoelastic  danmlng  In  an  Ideal  material.**  The 
equation  of  heat  conduction  In  his  case  Is  developed  by  a  slightly 
different  approach  from  that  used  by  Jeffreys  (6) .  Blot  considers 
the  unit  volume,  which  Is  fixed  in  space,  as  the  fundamental 
thermodynamic  system.  Unlike  Jeffreys,  however.  Blot  does  not 
find  It  necessary  to  use  the  conservation  of  mass  requirement. 

A  linearization  of  the  entropy  expression  leads  to  a  linear  heat 
conduction  equation  In  Blot's  treatment.  Quantities  of  transport 
by  the  transfer  of  mass  across  the  system  bottndary  are  not  con¬ 
sidered  In  (5)  and  do  not  arise  In  (6). 

In  this  report.  Intensive  quantities  per  unit  mass  are  used. 
The  boundary  of  the  fundamental  system,  the  unit  mass.  Is  varying, 
but  It  always  contains  the  same  mass.  The  usual  conservation  of 
mass  relationship  Is  used  and  the  equations  of  motion  are  shown  to 
be  consistent  with  the  unit  mass  system.  The  second  law  of  thermo 
dynamics  Is  postulated  In  deriving  the  rate  of  entropy  variation 
In  an  elastic  solid.  It  Is  then  shown  that  energy  dissipation  Is 
due  explicitly  to  the  Irreversibility  of  the  heat  conduction 

Srocess.  The  non-llnearlty  In  the  heat  equation  Is  then  shown  to 
e  necessary  In  order  that  predictions  of  the  theory  be  consistent 
with  the  energy  conservation  requirement.  This  Is  a  point  of  In¬ 
terest  In  that  a  linear  thermodynamic  consideration  leads  to  a 
non-linear  field  equation.  It  Is  also  a  point  of  some  importance 
In  the  calculation  of  energy  dissipation  per  cycle  In  the  steady- 
state  vibration  of  a  system  because  the  amount  of  heat  ejected  by 
the  system  Is  entirely  lost  by  the  linearization. 

An  example  ^Ich  Is  simple  enough  to  be  worked  out  In  mathe¬ 
matical  detail  Is  solved  and  the  energy  dissipation  calculated. 

An  approximate  theory  Is  then  Introdiiced.  The  res\slt  Is  shown 
to  be  very  accurate,  as  compared  to  that  obtained  by  the  exact 
method . 


*  Nuiid>ers  In  parentheses  refer  to  the  xsfsrmcM  on  page  30. 

**  Sometimes  the  treatment  of  thermoelastic  damping  In  ideal¬ 
ized  materials  is  erroneously  considered  as  to  explain  actual 
damping  in  real,  imperfect  materials. 

Manuscript  released  for  publication  1  April  1962  as  an 
ASD  Technical  Doeueentaxy  Rsport* 
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II.  THERMODYNAMICS  OF  ELASTIC  SOLIDS 


In  the  consideration  of  the  cottpled  phenomenon  of  mechanical 
deformations  and  the  conduction  of  heat  In  an  Idealized  elastic 
material,  ve  shall  choose  as  thermodynamic  state  variables  the 
six  components  of  the  stress  tensor,  <ru  ,  the  six  coiiq)onents 
of  the  strain  tensor,  <jj  ,  and  the  local  absolute  tenmerature, 
T.  We  are  dealing  with 'a  perfect  continuum  for  which  the  above 
state  variables  are.  In  general,  fimctlons  of  the  space  coor¬ 
dinates,  X{  ,  as  well  as  the  time,  t. 

An  Isotropic  solid  Is  defined  as  elastic  In  this  paper  If 
the  equations  of  state. 


o-j|  =  2/1  «||  +  X  -  (  3X  +  2/i.)  a  8jj  0  ,  (2-1) 

i  ,  j  ,  k  =  1 , 2 ,  or  3  , 

are  observed  everywhere  within  the  solid.  In  equation  (2-1),  6 
denotes  the  local  temperature  In  excess  of  the  vinlform  reference 
temperature  Tq  which  Is  chosen  together  with  the  reference  Tin- 
stressed  and  undeformed  state,  l.e.,  **1]  ~0^bTp;  X 

and  ft  denote  the  Lam6  constants  of  elasticity  under  Isotnermal 
conditions;  o  denotes  the  linear  coefficient  of  thermal  expan¬ 
sion  and 


8|j  =  I  if  1  =  )  , 

=  0  if  ijt)  . 

The  basic  relationships  to  be  studied  are  the  equations  which 
state  the  conservation  of  momentum,  mass,  and  energy,  and  the 
second  law  of  thermodynamics  for  the  Irreversible  process  written 
In  the  form  of  an  equality.  Intensive  quantities  are  to  be  used 
In  these  eqxiatlons.  These  quantities  are,  strictly  speaking, 
material  point  functions,  but  the  mlt  mass  may  be  used  to  help 
vlsTiallze  a  physical  thermodynamic  system,  e.g.,  the  value  of 
the  specific  entropy  at  a  point  may  be  regarded  as  the  total 
entropy  of  a  unit  mass  centered  on  the  point  In  question.  Since 
the  boundary  of  this  fundamental  system  Is  deforming  with  It, 
there  Is  no  mass  transfer  across  tne  boundary  and  the  question  of 
quantities  of  transport  need  not  arise. 

The  conservation  of  momentum  Is  stated  by  the  well-known 
equations  of  motion 
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+  /)  F,  «  /9  u,  . 


(2-2) 


"■jl.J 

In  the  last  equation  p  is  the  mass  density  of  the  solid;  Fi,  the 
bodily  force  per  unit  mass;  and  Uj[,  the  displacement  vector.  We 
have  also  adapted  the  convention  or  denoting  differentiation  with 
respect  to  space  coordinates  by  a  comma,  and  differentiation  with 
respect  to  time  by  a  dot  over  the  quantity.  It  Is  perhaps 
worthwhile  to  remark  that  equations  (2-2)  are  satisfied  by  the 
medium  In  Its  deformed  state,  and  hence,  the  use  of  Etilerlan 
coordinates  are  Indicated.  This  Is  consistent  with  the  use  of 
a  deforming  unit  mass  system.  In  addition  to  the  eqiiatlons  of 
motion,  appropriate  stress,  or  displacement,  boundary  conditions 
must  be  satisfied. 

The  conservation  of  mass  can  be  stated  In  the  simple  form 


where  c  Is  the  cubical  dilatation  and  Pq  Is  the  density  of  the 
material  at  zero  strain.  In  addition,  there  are  the  usual  com¬ 
patibility  requirements  on  the  strains  such  that  no  material  dis¬ 
continuity  may  occur. 

The  conservation  of  energy  may  be  stated  by  the  equation 


at 


(2-4) 


The  terms  In  the  above  equation  represent,  respectively,  from 
left  to  right,  the  rate  of  change  of  Interaal  energy,  the  rate 
of  change  of  specific  kinetic  energy,  the  rate  of  work  being 
done  by  the  bodily  forces,  the  rate  of  work  being  done  by  the 
surface  tractions,  and  the  divergence  of  a  heat  flux,  expressed 
for  the  volume  occupied  by  a  \inlt  mass  at  a  partlcvilar  time. 

The  expression  for  the  work  done  by  the  stirface  tractions  can  be 
derived  by  Integrating  the  scalar  product  of  the  sxirface  traction 
and  the  velocity  of  a  surface  element  over  the  boundary  of  a  unit 
mass,  and  by  applying  the  divergence  theorem  after  converting  the 
surface  traction  Into  stresses 
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/  T|  U|  </A  -  /  «r|j  u,  »,  </A 
•'A (unit  mast)  •'A 


</v 


We  can  eliminate  the  kinetic  energy  term  from  equation  (2-4)  by 
scalarly  multiplying  the  eqtiatlons  of  motion  (2-2)  by  the  vector 
U|  and  sxibstltutlng  the  results  Into  equation  (2-4),  obtaining 

Hi  =  ^  J_  j 

dt  p  b\  ~  p  '»•  (2-5) 

The  strain-defining  relationships 


) 


have  also  been  used.  The  substitutions 


(2-6) 


€ll  =  € 


1^22 


=  €, 


“^3  >2^23 


>2€3,-€5,2€|2*€6 


may  be  used  to  transform  eqiiatlon  (2-5)  Into  (using  differen¬ 
tials)  the  eqviatlon 


=  OTj  +  dq 


(2-7) 


with  d(l  Indicating  a  differential  heat  added  to  the  unit  mass, 
l.e.. 


(2-8) 
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The  equations  of  state  In  terms  of  the  new  variables  £,  and  Q* 
become  '  ' 


/»  oTj  «  2/x  C| -I- X  c  -  (  3X -f  2;a) a  d 
for  I  a  1 , 2  or  3  *, 


(2-9) 


and 


^o-,  for  i=  4,5  or  6  , 

(2-10) 

In  order  to  derive  the  expression  for  the  change  of  specific 
entropy,  denoted  by  the  symbol  s,  the  following  second  law  of 
thermodynamics  Is  postulated: 


Tdz^dq  ,  (2-11) 

the  relationship  for  a  reversible  process.  This  postulate  relies 
on  the  argument  that  the  only  type  of  "driving  force"  besides  the 
temperature  gradient  are  the  mecnanlcal  forces;  and  that  an  elastic 
body  eventutuly  comes  to  rest  under  a  sustained  load,  \inllke  the 
continuous  flow  of  heat  under  Its  "driving  force,"  the  temperatxire 
gradient.  The  behavior  of  an  elastic  body  Is  also  unlike  the 
viscous  flow  of  a  fluid  under  a  sustained  shearing  force. 

The  acceptance  of  this  relationship  as  applicable  outside  of 
thermodynamic  equilibrium  Is,  of  coxurse,  an  essential  part  of  the 
usual  theory  of  elastic  solids.  Also,  In  accepting  eqxiatlon 
(2-11),  we  assxime  that  the  specific  entropy  depends  explicitly  on 
the  state  variables  only. 

We  are  now  ready  to  derive  an  expression  for  the  specific 
entropy.  Writing 


flfq  Tj  +  c  dT 


(2-12) 


where  the  and  c  are  as  yet  undetermined  functions  of  the 

state  variables,  and  substituting  this  expression  Into  equations 
(2-7)  and  (2-11),  we  obtain 


</U  «(o-,  +  T,  )</€,+  Cfl^T  , 


(2-13) 
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and 


dJ 


(2-14) 


Since  U,  as  well  as  s.  Is  a  function  of  state  only,  dU,  as  well 
as  ds,  Is  a  perfect  differential.  This  will  furnish  sufficient 
relationships  to  determine  the  nature  of  Ti  and  c.  We  have  from 
equation  (2-13), 


d(a,  +  T.) 


and 


d(o"l  +  T|) 


^(O']  +  Tj) 

d«, 


ai 


_ac 

a«, 


and  from  equation  (2-14) , 

dr 

de 


j 


dr 

* 

d« 


and 


aj 


I  ac 

T“a7 


From  eqxiatlons  (2-16)  and  (2-18),  we  may  see  that 


dT  '  T  ' 


da 


1  . 


dT  dT 


(2-15) 


(2-16) 


(2-17) 


(2-18) 


(2-19) 


so  that 


T 


(2-20) 


which,  from  the  equations  of  state  (2-9,10),  Is  equivalent  to  say¬ 
ing  that 
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Ti  ®  ^(3X+2/t)oT  ,  for  i*l,2or3  (2-21) 

as  0  for  i  *  4, 5  or  6 

We  see  that  reqtilrements  (2-15)  and  (2-16)  are  satisfied  by  the 
T:  thus  determined.  In  addition,  from  equations  (2-16)  and 

(2-21), 

do  d  Ti 

l.e.,  c  Is  Independent  of  any  of  the  strains. 

The  differential  heat  added  to  a  unit  mass  can  itow  be  written 
In  terms  of  known  quantities 


/9flfq=(3X'i'2ft)aT^<  +  pc  dT 


(2-23) 


Hence  we  see  that  c  Is  just  the  specific  heat  of  the  material 
observed  under  constant  strain  conditions. 

By  assuming  Fourier's  Law  for  heat  condxictlon. 


J 


kT,, 


(2-24) 


the  field  equation  for  T 


kT,,, -/.e  (2-25) 

Is  finally  obtained.  Of  course,  k  Is  the  thermal  conductivity  of 
the  material.  The  Inclusion  of  the  last  term  differentiates 
equation  (2-25)  from  the  ordinary  heat  conduction  equation  de¬ 
rived  by  neglecting  the  thermal  expansion  of  the  material.  The 
non-llnearlty  of  this  extra  term  plays  an  Important  role,  as  we 
shall  see  later.  In  the  determination  of  thermoelastic  dissipation 
In  steady-state  vibrations  of  elastic  bodies.  We  shall  only  point 
out  here  that  a  non-llnear  field  eqxiatlon  has  been  derived  from 
linear  thermod;^amlc  considerations.  When  T  In  the  last  term  of 
equation  (2-25)  Is  taken  to  be  constant,  the  equation  forms  the 
basis  of  so-cailed  "linear  thermoelasticity." 
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The  heat  conduction  equation  (2-25),  the  emiationa  of  motion 
(2-2),  the  equations  of  state  (2-1),  t1:»  compatibility  require¬ 
ments,  and  the  appropriate  stress,  displacement,  and  temperature 
boundasry  conditions  form  the  basic  boundary  value  problem  of  the 
area  of  study  commonly  classified  under  the  heading  of  thermo- 
elastlclty. 

The  solution  of  this  problem.  In  the  linearized  version.  Is 
suggested  In  reference  (5),  along  with  the  now  well-known  varia¬ 
tional  principle  of  thermoelasticity  which  provides  a  powerful 
tool  for  the  numerical  evaluation  of  soltxtlons. 

We  note  that  for  a  given  problem  the  equations  of  motion  can 
be  solved.  In  principle  at  least.  In  terms  of  the  temperatures. 
The  dilatation  <  can  then  be  substituted  Into  the  heat  conduc¬ 
tion  eqtiatlon  which  will  contain  the  tenperature  as  the  only 
unknown. 

Pl^slcally,  we  see  that  a  tenmerature  variation  gives  rise 
to  a  dilatation  which  plays  the  role  of  a  heat  source,  causing  a 
change  of  the  temperat\ire  and  fxnrther  deformations. 


III.  THERMOELASTIC  DAMPING 


Inasmuch  as  all  thermoelastic  processes  Involve  the  conduc¬ 
tion  of  heat,  they  are  dissipative  In  nature.  We  shall  discuss 
this  point  In  detail. 


38  Let  us  consider  the  rate  of  change  of  specific  entropy 

.  Equations  (2-8)  and  (2-11)  can  be  combined  Into  the 

equation 


P 


^  — -  J_ 

at  T 


j 


which  may  be  written  in  the  form  of  a  so-called  "balance 
equation:"  (7) 


(3-1) 


as  _  ,  I 

'’aT  ‘t 


H.' 


(3-2) 


The  chaxige  of  the  specific  entropy  Is  seen  tc^come  from  two  causes: 

the  negative  divergence  of  an  entropy  flow  -iji^  and  an  entropy 

production  whose  strength  Is  given  by  the  last  term  of  the  balance 

dim  —.1  ^8t 


equation.  Denote  these  two  parts  by  _ _ 2L  and  « 

^  d\ 


d\ 


8 


respectively,  then 


dt  at  at 


(3-3) 


The  local  entropy  production  Is,  In  the  sense  o£  general 
Irreversible  thermodynamics,  the  prod^uct  of  a  "driving  force," 


and  a  flux,  In  this  case  the  heat  flux  J. 
phenomenological  lav  Is  assumed,  l.e.. 


(3-4) 


If  a  linear 


J,  =LX, 


(3-5) 


we  obtain  Fourier's  Law  of  heat  conduction,  usually  written  In 
the  form 


J,  =  -kT,,  .  (3-6) 

L 

The  conductivity  k  Is,  therefore.  Identified  with  ~  jt 

Equation  (3-6)  has  been  used  In  the  derivation  of  the  teiiq>erature 
field  equation  (2-25). 

The  local,  or  Internal  entropy  production  Is  a  measure  of 
the  Irreversibility  of  the  process.  It  Is  a  positive  definite 
quantity.  We  see  from  equation  (3-2)  that  the  production  of 
entropy  In  the  case  of  thermoelastic Ijtv  originates  from  the 
presence  of  a  temperature  gradient.  The  cause,  and  the  only 
explicit  cause,  of  dissipation  Is,  therefore,  the  flow  cf  heat 
from  a  higher  tenq>er attire  toward  a  lower  tenqierature. 

In  a  stationary  state  (7),  the  entropy,  being  a  function  of 
state,  remains  constant.  The  first  member  of  equation  (3-3) 
vanishes.  Since  the  same  amount  of  heat  Is  entering  the  system 
at  a  higher  temperature  as  la  leaving  at  a  lower  tenqierature, 
more  entropv  Is  flowing  out  of  the  system  and  this  amount  must 
be  made  up  locally  In  order  to  fulfill  equation  (3-3). 
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We  shall  direct  our  attention  to  a  specific  aspect  of  the 
general  analysis  of  Irreversible  processes  Involving  elastic 
solids.  This  Is  the  area  which  may  be  loosely  described  by  the 
term  thermoelastic  damping.  To  dlstlnwlsh  It  from  the  study  of 
daiiq>lng  In  real  materials,  we  remark  that  thermoelastic  damping 
does  not  Inclvide  the  study  of  the  many  known  causes  or  mechanisms 
of  "conversion"  of  work  Into  heat  flow,  for  which  bridged  micro¬ 
scopic  and  macroscopic  theories  are  needed.  Rather,  It  comprises 
a  study  of  what  happens  dixrlng  heat  flow  due  to  mechanical  de¬ 
formations  In  an  Idealized,  perfect  continuum.  Looking  at  It 
from  a  different  viewpoint,  this  study  Is  directed  toward  an 
answer  to  the  question:  much  of  the  total  danq>lng  capacity 

Is  due  to  microscopic  imperfections  In  a  given,  real  material? 

To  be  more  specific,  let  the  volume  occupied  by  an  elastic 
body,  and  bounded  by  the  s\irface  S,  be  V.  Let  this  body  be 
Initially  In  thermodynamic  equilibrium,  l.e., 

<^1J  =  «I1  -0  , 

and  T  =  Tg 


at  time  t  =  0 


For  t  >  0,  periodic  bodily  forces  and  surface  tractions  are  ap- 

?lled.  If  the  boundary  conditions  for  the  temperature  are  tlme- 
ndependent,  and  If  steady- state  response  can  eventtially  be 
reacned,  we  shall  attempt  to  determine  the  net  amount  of  mechan¬ 
ical  energy  required  for  each  cycle  to  maintain  this  steady-state. 

The  system  defined  above  Is  analogous  to  a  "heat  pump." 

The  working  medium,  the  elastic  solid.  Is  adjoined  to  an  Infinite 
reservoir  at  constant  teimeratvure  To  through  condtictlng  and 
Insulating  boundaries.  When  this  system  Is  being  "operated"  by 
the  action  of  the  applied  forces,  a  tenqperatxire  field  different 
from  To  will  be  observed  In  V,  Inducing  heat  flow  and  Increasing 
the  total  entropy  of  the  system  plus  the  resein^olr. 

Hence,  any  operation  of  this  pimtp,  excepting  pure  shear, 
will  be  an  Irreversible  process.  A  net  amount  oi  positive  work 
must  be  done  on  the  system  even  In  steady-state  operations  where 
the  temperature  and  the  strains  return  to  their  respective 
values  after  each  cycle.  From  the  first  law,  we  may  conclude  that 
a  net  amount  of  heat  must  be  "pumped"  out  of  the  system  for  each 
cycle.  Therefore,  we  can  Immediately  deduce  that  that  the  average 
tenq>erature  over  a  cycle  must  be  raised  above  T^  .  This  Is 
accoiiq>llshed,  of  course,  during  the  transient  stages  before  the 
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steady-state  Is  achieved. 

The  work  done  by  the  applied  body  forces  and  surface  trac¬ 
tions  per  cycle  Is  given  by 


-(3X  +  2M)a8|l]</v</t  /■^lpu,}aycf1 


Also 


T|  =  », 


on  s 


(3-8) 


In  the  above  eqixatlons  ot  is  the  angular  frequency  of  excita¬ 
tion  and  Is  the  unit  normal  vector  on  S,  pointed  outward. 
By  the  divergence  theorem. 


i/, 


dM 


dM 


and  by  the  equations  of  motion  the  last  term  may  be  combined  with 
the  Integral  of  the  work  done  by  the  body  forces. 


Fj  U|  dy 


U|  dy 


4^ 


ip  ii^)dy 


(3-10) 
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The  equations  of  state  have  also  been  used  In  deriving  equation 
(3-7). 


Assuming  the  deformations  are  small,  V  and  S  may  be  consid¬ 
ered  Independent  of  time,  and  the  order  of  Integration  may  be 
Interchanged  In  equation  (3-7), 


€„)  + 


_a, 

at 


(€„€ 


dy 


at 


( p  u  I )  dy 


(3-11) 


Because  of  periodicity,  the  first  and  the  last  Integrals  vanish 
and 


w«==-^  J^(3\-h2/JL)a6  j^dy  di  . 


(3-12) 


The  total  heat  flow  out  of  V  per  cycle,  Hq,  can  be  evalu¬ 
ated  by  keeping  accoimt  of  the  flow  across  S.  It  Is  given  by  the 
expression 


Ho=-k/  /  6,,  V,  di  d\  . 

M  •'s 

Using  the  divergence  theorem  «^aln,  we  obtain 

From  eqxiatlons  (^2-25)  and  (3-14), 

,t  + 


(3-13) 


(3-14) 


X  [‘3X.2^)a(T..0)  -|i *p C |f]avrft  0-15) 
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or 


Ho  **  J^(dX+ 2^)  ad dv  c/t  .  (3-16) 

Thus  we  have  the  Identity 

w,  =  H,  (3-17) 

as  expected. 

An  Interesting  point  has  been  brought  otit  In  the  above 
derlvat'lon.  We  note  that  If  eqtiatlon  (z-25)»  the  heat  conduction 
relationship  had  been  linearized,  l.e..  If 

kT,',,  c  |^  +  (3X  +  2/LL)aTo  (3-18) 

Is  used  In  equation  (3-15),  we  obtain.  Instead  of  eqtiatlon  (3-16), 


(3-19) 


The  value  of  Wq  Is  not  affected  by  such  linearization.  We  are  led 
to  an  apparent  contradiction  to  tne  first  law. 


W.ikH;  .  (3-20) 


In  equations  (3-18),  (3-19),  and  (3-20)  the  prime  Is  used  to 
associate  the  quantity  In  question  to  T'.  the  teoqieratxire  field 
defined  by  the  linearized  equation  (3-18).  Since  the  absolute 
temperature  T  Is  not  expected  to  vary  appreciably  from  Tp,  (as 
the  thermal  boundary  conditions  are  time  Independent),  the  non¬ 
linearity  In  equation  (2-25)  Is  Indeed  small,  and  T'  shoxild  be 
a  good  approximation  for  T.  But  as  Just  point  out,  one  must 
exercise  care  If  the  energy  dissipation  per  cycle  of  steady-state 
vibration  Is  sought. 
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IV.  AN  EXAMPLE 


The  exact  evaluation  of  thermoelastlc  energy  dissipation  re> 
quires  solving  simultaneously  for  the  stress  and  tenq>erature 
fields  from  the  coupled  equations  of  motion  and  heat  condtictlon. 
We  shall  demonstrate  some  of  the  typical  procedures  Involved  In 
such  an  evaluation  ^^th  the  aid  of  a  one -dimensional  exaitqile  of 
an  Infinite  plate  with  uniform  normal  pressures  applied  sinus¬ 
oidally  on  both  flat  surfaces.  From  the  solution  of  this  example 
some  lnq>ortant  principles  are  dedticed,  leading  to  the  development 
of  a  general  theory  In  Section  V. 


Let  us  choose  the  coordinate  system  such  that  the  plane 
X|— 0  coincides  with  the  mid-plane  of  the  plate.  If  the  plate 
thl'ckness  Is  2h  and  the  applied  pressure  Is  P  sin  u  t,  the  plate 
will  be  deforming  In  an  extension-contraction  mode  perpendicular 
to  the  plate.  The  lowest  natural  frequency  of  the  plate  vibrating 
In  this  mode,  determined  by  assuming  Isothermal  conditions.  Is 
given  by  the  well-knowii  expression 


ir  / 

“  2h  V  P 


(4-1) 


In  the  example,  we  will  assume  that  the  temperatiires  on  both 
surfaces  of  the  plate  are  to  be  kept  at  the  constant*  reference 
temperature  T,  for  all  time.  Since  all  the  quantities  are  func¬ 
tions  of  time  md  X I  only,  and  since 


U 


U. 


=  0  , 


(4-2) 


the  strains 


^22*^33*^23=^31*^12=0  * 

and  £  ,1 

^11  • 

The  cubical  dilatation  Is,  therefore,  simply 


(4-3) 

(4-4) 


The  stresses  are 


(4-5) 
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^*12  “^23“®’3I  > 

(4-6) 

0’22®®’33®^*  “(3X+2/A)ad, 

(4-7) 

and 

<r||  »(X*f2fL)< -(3X-f2/i)ad . 

(4-8) 

The  equations  of  niotlon(2-2)  are  reduced  to  only  one, 


(4-9) 


In  tdilch  the  variation  of  the  density  Is  assisned  small.  By 
limiting  ourselves  to  the  case  where  the  exciting  frequency  oi  Is 
well  under  (Dq  ,  say  ,  the  right-hand  side  of  equation 

(4-9)  Is  extremely  small,  laqplylng  that  Is  almost  Inoe- 

pendent  of  X|  .  Let  us  assume  that 


=  0 


Since  0*11  Is  eqtial  to  the  applied  loads  on  the  boundary.  It 
must  be  so  everywhere  within  the  plate  thickness. 


<r|  I  =  P  sin  a»t 


(4-10) 


This  means,  from  equation  (4-8),  that 

P 


<  — 


sinwt  +  ^  a  6 

X+2/1  X  +  2/i 


The  heat  conduction  equation  (2-25)  can  now  be  reduced  to 


k  A  -/.c  —  +(3X+2,i)o(  V 
dx  d  t 

X  +  2/i  «t  J  _ 


j  Poi  COSwt 

L  X+  Zfi 


(4-11) 


(4-12) 
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The  solution  6  must  also  satisfy  the  botmdary  conditions 


0(±h,t)  =  O 


(4-13) 


for  all  t.  Since  we  are  seeking  for  the  steady- state  solution  of 
6  ,  the  condition  of  periodicity 

0(X,  ,t)*d(X,,t  +  2ir/u,)  (4-14) 

must  be  applied  for  all  values  of  X  |  within  the  plate. 

We  may  non-dimensional ize  eqtiatlon  (4-12)  by  Introducing  the 
dimensionless  variables  C  >  t  and  ^  which  are  defined  as 
follows : 


^  h 

(4-15) 

T  —  Cl>  t  , 

(4-16) 

and 

. 

■o 

(4-17) 

Upon  substitution,  we  obtain  the  equation 

d<b  dit  dS 

dc  at  •  '  dr 


(4-18) 


The  coefficients  C|  >  0]  »  and  b|  y  are  also  non-dimensional  and 


c,  * 


p  c  hw 


and 


^  ^  P(3X-»-2M)ah<a 
'  (X  +  2/i)k 

(3X-»-2Mi*a*Tnh*a» 


‘  (X  +  2/i)k 

The  boundary  condition  (4-13)  In  terms  of  the  new  quantities  Is 


(4-19) 

(4-20) 

(4-21) 
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^(±l,T)-0  , 

and  the  periodicity  condition  (4-14)  Is  now 


(4-22) 


^(C,T)*^(C,r  +  2ir)  (4-23) 


for  all  C  and  r 

For  an  Idealized  elastic  material  having  physical  properties 
similar  to  the  common  carbon  steels,  the  ratios  among  C|  >  Qi  > 
and  b|  asstme  the  proportionalities 


c,:  a,:  b|»  400  :l  :l 


(4-24) 


when  the  presstxce  anmllttide,  P,  Is  as  high  as  14,000  psl,  actually 
equivalent  to  more  than  a  half  of  the  stress  mamltude  needed  to 
cause  the  material  to  yield.  Hotlce  that  equation  (4-24)  holds 
Independently  of  the  quantities  k,  h,  and  .  Since  the  Induced 
temperature  change,  0  ,  Is  not  expected  to  be  large.  It  Is 

convenient  to  miuce  the  transformation 


(4-25) 


In  eqtiatlon  (4-18)  for  the  purpose  of  comparing  relative  magnitudes 
of  tM  terms  on  t^  right-hand  side  of  the  equation.  The  magnitude 
of  ^  Is  now  of  the  order  tmlty.  Then 


d?^*®'“d^ 


d4r 

+  COSr  +a,(  I  +/9^)  — 


where  /3 


(4-26) 

(4-27) 


has  been  used.  The  boundary  and  periodicity  conditions  are  now, 
respectively. 


^(±I,t)  «  0  , 


(4-28) 


17 


and  +2ir)  . 


(4-29) 


Since 

l^'l  «  • 

all  non-linear  terms  in  e<tuation  (4-26)  may  be  considered  small. 

The  non-linear  botmdary  valvie  problem  defined  by  equations 
(4-26),  (4-28),  and  (4-29)  will  be  solved  by  a  perturbation 
scheme.  The  method  id.ll  depend  on  the  smallness  of  the  two 
_  a, 

quantities  p  and  -jf-  (as  coiiq>ared  to  vmity).  Let  us  expand^ 

C| 

into  a  power  series  of  the  parameter  p  ; 

=  I  \KC,t)/9"  (4-30) 


and  then  substitute  into  both  sides  of  the  differential  equation 
(4-26)  to  obtain 


. ) 

~  t  . COST 

+  + . ) 

(4-31) 

+  C|  /9  (  ^0  . ^  COST 

+  . )  ^  (  ^  ^2  . ^  • 


Terms  containing  like  powers  of  /3  on  both  sides  of  the  last 
equation  may  be  collected  and  equated,  giving  a  set  of  linear 
differential  equations  for  the  coefficients  of  the  power  series. 


. Y  *  (  Cl  +0|  )  ■*'  Cl  COST 

'  '  dT  ' 


(4-32) 
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(4-33) 


dV 


d^i 

(c,+a,)—  +  c,^(r^co8T  +a,^o”J7 


I 


(c, 


d% 

+a,  )-^  +  c,  COST  +  a, 

dr 


9 


(4-34) 


The  problem  Is  therefore  reduced  to  solviim  the  above  set  of  dif¬ 
ferential  equations  in  the  given  order.  The  boundary  and  perio¬ 
dicity  conditions  are,  respectively. 


^„(±I,t)*0  (4-35) 

and  (€,  f)  =  ,  t  +  2v)  (4-36) 

for  each  of  the  coefficients. 

For  n  •  0,  the  solution  may  be  found  by  writing 

—  f^(^)  COST  +  Sinr  ,  (4-37) 

where 

c'l  C,  +  0,  .  (4-38) 

and  substituting  Into  equation  (4-32).  Equating  terms  which 
contain  the  common  factor  cos  t  ,  we  get 

«oK)-c'|9.(€)=0  ; 


equating  terms  i^lch  contain  the  common  factor  slnT  we  get 

g;({)  +  c;f,(€)-0  ,  (4-40) 
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The  boundary  conditions  are 


'a(±l) 


9„(±l) 


The  solutions  for  and  g  are 

o  o 


0 
S' 


(4-41) 

(4-42) 


fo(^ )  =  Ao  siny C  sinh^/ -S-C  -  Bq  cos  ^  cosh^  ^ 


(4-43) 


and 


with 


C  +  B^sin, 


€  sinK 


e  ,(A-44) 


and 


(4-45) 


(4-46) 


which  may  be  verified  by  direct  siibstltution. 

In  the  case  of  low  carbon  steel,  we  find  for  the  non-dimen¬ 
sional  coefficient  C|  the  value  of  4.5  x  10‘,  taking  h  >  1  inch 
and  u  "  30X  dtp  .  Based  on  this  knowledge,  we  can  make  the 
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following  estimate  of  the  magnltxide  of  f^  (  C  ) 


and,  similarly,  for  C  ) 

(4-48) 


S  to(f )  +  g*(€)  +  2  4  I  «o(<)|+ 14  f=4+4+4+|.l3,(*-«) 

c,  c, 

or  \%\  ^3.36  approximately. 

The  magnitude  of 


is,  therefore,  limited  to 

13.36 X/8|  i  0.0084  , 
or  9  S  0.0084X273**  2.2®  F  . 

Qi  I 

Since  the  fraction  ^  is  about  ,  and  with  the  above 

limitations  on  the  ma^ittide  of  %  ,  the  differential  e<iuatlon 

for  ^1  may  be  simplified  by  omitting  terms  which  are  proportional 
to  Qi  when  the  same  functions  appear  with  Ci  as  a  factor. 
Remembering  that  the  solution  for  will  be  multiplied  by  the 
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factor  /9  in  the  final  aolutlon  tf/  ,  It  la  reasonable  to  as- 
aune  that  this  qn>i^oxlmatlon  irlll  not  Introduce  serious  errors  In 
the  total  solution.  With  this  simplification,  equation  (4-33) 
becomes 


.  if, 

__C,— +c,f„co.T 


(4-51) 


The  solution  Is 


% 

vhere 


+  ^i(0  co8  2t+ 

C|  L  4c,  . 


8in2T  (4-52) 


f,  (C)  *  A,  8ln^/c[C  8inh-yc[ €  -  B,  co8>/c}  C  C08h/c5  €  j  (4-53) 


Q,(C)*  A,  008-/?, C  C08h-w^^  +  B,  8ln-/^ ^  slnh-/c[C  ,  (4-54) 

with 


Cf  008/^  C08h-/H 

A,  = - ; -  ;  (4-55) 

4c',  ( 8inhVcJ  +  co8Vc[) 

c*,  8in-/c5  8lnhyc[ 

B, « - - - ! -  .  (4-56) 

4  c',  ( 8inh*/c[  +co8Vc^  ) 

Additional  perturbations  may  be  obtained  in  like  manner.  They  are 
not  necessary,  however,  due  to  the  rapid  convergence  of  the  series 
of  powers  of  p  .  The  mamltudes  of  f ,  ,  Q ,  ,  and  vi  can  be 
estimated  and  are  comparable  to  ^  and  g  .  We  may  summarize  by 
stating  that  the  aolutlon  Is  very  accurately  represented  by  taking 
onlv  the  leading  term  of  the  aeries  (4-30),  and  that  there  appears 
to  be  no  need  to  go  beyond  a  two-term  approximation  for  the 
solution. 
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The  HDXk  done  during  a  cycle  la  given  by  the  general  expres¬ 
sion  derived  in  Section  III: 


Wo«-(3X  +  2/i)a 


♦  4.  SI 

■I  f 

•/f  Jy 


d^cfy  d\ 

at 


(4-57) 


Choosing  the  voltaae  with  a  unit  area  perpendicular  to  the  X|  -axis 
as  the  volume  V,  and  Inserting  expression  (4-11)  for  «  into  the 
above,  we  have 

^  _  (3X-l‘2jx)a P«tf 

®  ~  X  +  2fi  /  /  ® 

•'t  •'-h 


(3X-t-2/i)VX)  P*h 
(X  +  2/ifpc 


I 

Cfr 


(4-58) 


It  is  seen  that  only  those  terms  uhich  contain  the  factor  cost 
in  the  solution  for  ^  contribute  toward  damping,  excliiding  all 
terms  except  ia(  C  )  t  energy  dissipation  is, 

therefore,  given^y 


(3X-t-2/tfa*ToP*h 
(X  +  2/i)/)C 


I 

fo(0</C  dr 


(3XH-2ft)*a*ToP*hv  _Cl  (sinhyS?,  -sin^/gc;) _ 

(3X  +  2/t)‘o*ToP*hir 

(X+2/t)Vcy^ 
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As  a  check,  fornula  (3-14)  may  be  used  to  compute  the  heat 
flow  across  the  two  unit  surface  areas, 

.t+** 


V 

T+2ir  m\ 


dr 


-nT.g  r  f  ^ 

h«  Jr  J-i  4{* 

By  eq^tlon  (4-30), 

. )« 

Using  equations  (4-32),  (4-33),  (4-34),  etc.. 


dr 


,T+2ir  ^1 


r  -I 

kTJS  r  r  dfc 

H„ - jj— |^(c,  +  a,)— +c,C0ST  J  a(clr 


kT^r^r'r  a*, 


COST+oUr 

•  0  d 


dr  . 


+  0.) 


^^2 

dr 


] 


(4-60) 


Since  the  functions  ,  4^1  ....  are  all  periodic,  we  see  that 

will  be  sero  \inless  at  least  a  two -term  approximation  for  if 
is  used.  In  ehich  case,  the  net  heat  flow  out  would  be  exactly 
equal  to  the  work  done  in  a  cycle.  To  arrive  at  a  basis  of  com¬ 
parison  with  the  actual  danq>lng  capacities  of  a  real  material  In 
this  mode  of  deformation,  we  express  as  a  percentage  of  the 
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maximum  stored  strain  energy,  E,  which  Is  reached  during  a 
cycle 


E 


‘ij 


</X, 


X+2/a 

hence 

(3X-»-2ft)VTe  ir 

ik+zfi)  pc,y^ 


(4-61) 


(4-62) 


which  Is  approximately  0.02Z  for  h  -  1"  and  w  *  30%  Wo 
This  Is  less  than  the  damping  which  wotild  occur  In  a  real,  liq>er- 
fect  material.  An  Interesting  problem  would  be  to  formulate  a 
phenomenological  law  of  Internal  danq>lng  In  terms  of  the  entropy 
production,  which  would  probably  lead  to  a  heat  conduction  eqxia- 
tion  with  additional  source  terms.  The  solution  might  not  be  any 
more  difficult  If  the  same  kind  of  approximation  ^ich  Is  used 
In  our  example  Is  made.  Care  must,  of  course,  be  exercised  so  as 
to  maintain  the  same  stress- strain-temper  attire  relationship. 


Let  us  return  to  the  expression  for  thermoelastic  daiqilng, 
equation  (4-59).  We  note  that  Is  dependent  of  ^e  frequency 
or  excitation  only  through  the  parameters  C|  and  C|  ,  and  that 


Mm  Wo  ~  M m  W©  »  0  (4-63) 

w—o  <*»—00 


Wo  is  a  positive,  continuous  function  of  the  exciting  fre¬ 
quency.  This  means  that  there  exists  a  frequency  at  which  the 
energy  dissipated  per  cycle  Is  at  a  maximum.  The  factor 


Is  shown  in  Figure  1  on  page  31.  We  see  that  the  maxinnm  is  reached 
at  a  very  low  frequency,  0.01  for  this  example,  and 

that  It  Is  reached  within  a  very  nainrow  frequency  band. 

V.  GENERAL  THEORY 


Within  the  range  of  linear  thermodynamics,  the  dilatation  < 
may  be  written  as  the  sum  of  two  parts:  the  fractional  change  of 
volume  due  to  the  temperature  and  the  fractional  change  of  voltime 
due  to  the  stresses,  l.e.. 


€«€'(d)  +  <" 


It  Is  known  that  s' (  0  )  Is  linearly  related  to  0 


(5-1) 


.'(6)  =  >:  8  +  9  .  (5.2 

where  0  Is  the  average  tenq>erature  over  the  entire  volume  at 
any  given  Instant;  and  X'  <ind  u  are  constants  containing  the 
elastic  modiill. 

The  general  heat  equation  (2-25)  may  now  be  written  In  terms 
of  «'  and  «"  : 


dt 


+  (3X  +  2,»).T,(|i  + 


8  8«“,  a.' 

To  8t  df 


To  dt  (5-3) 


We  shall  regard  all  terms  In  the  parentheses  as  heat  sources. 

The  solution  0  Is  known  to  be  small.  Let  us  examine  the  equa¬ 
tion: 


k7*0=/>c-|^+(3X+2/i)aTo 


To  dt  dt  TodtJ(5-4) 


In  which  0  denotes  a  known  function  of  small  magnitude  when  com¬ 
pared  to  To  •  The  steady- state  solution,  0  ,  of  equation  (5-4) 

will  be  mainly  due  to  the  first  of  the  heat  sources, 

(3X  +  2/»)aTo4^  . 

d  t 


By  neglecting  all  other  sources,  we  obtain  an  approximate  heat 
equation 
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de 


i£ 


It at 

^Ich  Is  linear,  and  which  Is  uncoupled  from  the  equations  of 
motion. 


(5-5) 


The  energy  dissipation  per  cycle  Is  a  steady-state  condition 

» 


The  energy  dissipation  per  cycle  is  a  steady- 
Is  given  by  the  double  Integral  [equation  (3-12)] 

►  t  +  Zw/<a 


Wo=(3X+2/i)a 

»  (3X+2/i)a 

since 

,  t  +2vAi» 


MTTCWfW  A 

J,  Vir 


e/yi  d\ 


(5-6) 
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)  d>/  d\ 
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Q  —  Q  from  periodicity  considerations 

at 

,t+2v/w 


W.  =-(3X  +  2/i)a 


/D  •  A 

,/  /e^rfvd. 

•/f  o/w  at 


(5-7) 


Multiply  both  sides  of  equation  <3-5)  by^/T^  and  Integrate 
over  V  and  hr  !  to  .  We  see  that  the  last  Integral  becomes 

^t+2v/w  -  -t+2v/w  - 

(3x+2/i)a  /  /  e  It  <'*  =  4- /  he„,a^at. 

•'t  Vy  Tj«'t  •'V  (5-8) 

By  the  use  of  the  divergence  theorem, 

f  ee.^fdy  ds^f  (9, 

Jy  •'s  •'V 


(5-9) 


Equation  (5-8)  Is  transformed  Into  an  approximate  formula  for  the 
quantity  of  work  done  per  cycle  of  steady-state  deformation: 
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(9,,)Vv</t  . 


(S-10) 


Since  either  B  or  ^,1  vanishes  on  the  boundary  S  by  the 

stationary  boundary  condition  assumption. 


By  means  of  equation  (5>10),  the  computation  of  thermo* 
elastic  damping  In  the  general  case  Is  reduced  to  a)  finding  the 
dilatation  by  assuming  Isothermal  conditions,  b)  determining  the 
temperature  field  from  equation  (5-5),  and  cl  evaluating  the  work 
done  bv  formula  (5-10).  This  Is  a  considerable  sliimllflcatlon 
over  the  "exact  method"  which  has  been  Illustrated  by  the  exasmle 
In  Section  IV.  In  that  example.  It  may  be  shown  easily  that  the 
above  approximate  formula  will  lead  to  a  r^stilt  which  la  the 
same  as  the  "exact"  formula  except  that  c.  Is  replaced  by  C| 

In  ecniatlon  (4-59).  It  may,  therefore,  be 'said  that  equation 
(5-10)  furnishes  a  much  sliq>llfled  route  to  calculatltig  the 
amount  of  energy  dissipation  for  small  deformations. 


VI.  CONCLUSIONS 


The  following  conclusions  are  reduced  from  observing  the  re¬ 
sults  and  their  derivations  In  the  above  discussion: 

(1)  The  field  equations  of  thermoelasticity  are  derived  by 
considering  the  conservation  of  mass,  momentum,  and 
energy  In  addition,  a  second  law  postulated  for  perfect 
elastic  materials  Is  found  necessary.  These  equations 
represent  a  phencmenologlcal  theory  consistent  with  the 
notions  of  continuum  mechanics  In  the  study  of  elastic 
deformations  and  heat  flow. 

(2)  Thermoelastic  damping  In  Idealized  materials  Is  the 
explicit  resiilt  of  the  Irreversible  character  of  heat 
conduction  only. 

(3)  For  real  materials.  It  Is  suggested  that  the  numerous 
macroscoplcally  observed  damping  laws  may  be  represented 
by  source  terms  In  the  heat  condiiction  equation. 

(4)  The  heat  conduction  (■'•uation  for  perfect  materials  Is 
ixni-llnear  althoiigh  a  linear  thermodynamic  process  has 
been  considered. 
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(5)  An  inconslstnncy  will  apMar  between  the  eolutione  of  the 
field  equations  and  the  first  law  if  the  heat  equation  is 
linearised  in  the  nanner  o  f  the  usual  themoelasticity 
theories  in  the  literature. 

(6)  In  processes  where  thensal  boundary  conditions  are  sta¬ 
tionary  in  time,  the  general  expression  for  the  ener^ 
dissipation  for  a  steady- state  cycle  of  deformation  is 
derived. 

(7)  A  sinple  one-dimensional  exaomle  is  given  in  Section  IV. 
The  equations  of  motion  (equilibrium)  are  solved  in  terms 
of  the  temperature  field.  iMch  is  in  turn  solved  from  the 
non-linear  heat  conduction  equation  by  a  perturbation 
scheme.  For  an  idealised  elastic  material  having  prop¬ 
erties  similar  to  low  carbon  steel,  the  said  non-lineiurity 
is  shown  to  be  unnecessary  only  if  the  energy  dissipation 
is  computed  by  the  work  done  on  the  system. 

(8)  In  the  given  example,  the  Induced  tenq>erature  variation 
is  small  (less  thim  2.20F  for  a  maxtanm  stress  of  14,000 
psi,  and  at  a  frequency  which  is  30  percent  of  the  first 
natural  freouency  of  the  system),  lue  amount  of  energy 
dissipated  during  a  cycle  is  also  extremely  small  (alx>ut 
0.02  percent  of  the  maxlmmi  strain  energy  reached  during 
a  cycle).  Based  on  these  facts,  the  heat  equation  is 
linearized  and  uncotq>led  by  neglecting  all  source  terms 
except  that  representing  tra  rate  of  change  of  dilatation, 
calculated  by  assuming  isothermal  deformations. 

(9)  The  above  sinpllfication  is  generalized  to  derive  an 
expression  for  the  energy  dissipation  in  arbitrary  bodies 
under  elastic  deformations. 
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